Bogoliubov Theory of acoustic Hawking radiation in Bose-Einstein Condensates 
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We apply the microscopic Bogoliubov theory of dilute Bose-Einstein condensates to analyze quan- 
tum and thermal fluctuations in a flowing atomic condensate in the presence of a sonic horizon. For 
the simplest case of a step-like horizon, closed-form analytical expressions are found for the spectral 
distribution of the analog Hawking radiation and for the density correlation function. The peculiar 
long-distance density correlations that appear as a consequence of the Hawking emission features 
turns out to be reinforced by a finite initial temperature of the condensate. The analytical results 
are in good quantitative agreement with first principle numerical calculations. 
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I. INTRODUCTION 

Thanks to the impressive advances in the cooling and 
manipulation of ultracold atomic gases, experiments are 
now able to address novel regimes where the physics of 
coherent matter wave is strongly affected by zero point 
quantum fluctuations [2, 0, 0, El 0] • I n particular Bose- 
Einstein condensed (BEC) atomic vapors appear as a ver- 
satile and efficient tool for observing a very fascinating 
manifestation of quantum fluctuations, namely Hawking 
radiation from acoustic black holes, the so-called dumb 
holes 0,0. 

Hawking radiation is a most celebrated, yet still un- 
observed prediction of quantum field theory on curved 
space-times, which consists of the conversion of vacuum 
fluctuations into observable radiation in the vicinity of a 
black-hole horizon H . Elaborating on the mathematical 
analogy between the propagation of sound waves in in- 
homogeneous and moving media and the propagation of 
quantum fields on a curved spacetime background, Unruh 
predicted in 1981 the occurrence of an analog Hawking 
emission of sound in any system showing a sonic horizon, 
i.e., an interface between a sub-sonic and a super-sonic 
region [§]. 

Very recently, the experimental realization of a dumb 
hole-like configuration in a flowingatomic Bose-Einstein 
condensate was presented in Ref.[14|. From a more gen- 
eral standpoint, a dumb- hole configuration for surface 
waves on a tank of moving water was realized in [lol ]: 
the reported observation of the conversion of positive- 
frequency waves into negative frequency ones can be con- 
sidered as a classical analog of the Hawking effect [IT| . 
The observation of a horizon in a microstructured optical 
fiber was reported in [12] ■ The possibility of simulating 
acoustic black holes using a ring-shaped chain of trapped 
ions was considered in [1 31 ] - 

So far, much of the theoretical work on the analog 
Hawking radiation in atomic Bose-Einstein condensates 
has used some gravitational analogy to characterize the 
properties of the emission. Our point of view is differ- 
ent and aims at developing a microscopic understanding 



of analog Hawking radiation starting from the general 
theory of quantum fluctuations in condensed matter sys- 
tems and without referring to any gravitational analogy. 
A similar approach has been adopted in the few last years 
by other researchers: A first step in this direction can be 
found in Ref. Q where a general theory of analog Hawk- 
ing radiation based on Bogoliubov theory is presented 
but no detailed analysis of the observable consequences 
is carried out. This point of view has been pushed fur- 
ther in the very recent papers [13, EH , where the general 
framework is developed in full detail and numerical re- 
sults for some observable quantities are also discussed for 
the case of a smooth sub- to super-sonic flow transition. 

In the present paper we report a completely analytical 
study of Hawking radiation from dumb-holes in atomic 
Bose-Einstein condensates. Our calculations are based 
on a direct application of the standard Bogoliubov the- 
ory of dilute condensates and no explicit reference is ever 
made to the gravitational analogy. In order to make 
the problem analytically tractable, we consider a sim- 
plest step-like configuration where the transition between 
the sub-sonic to the super-sonic regions occurs on a very 
short length scale. Even if the surface gravity of this con- 
figuration is formally infinite, still a thermal-like Hawk- 
ing emission is found at a temperature fixed by the heal- 
ing length. Closed analytical formulas for the density 
correlations are extracted, which are found in excellent 
agreement with the numerical simulations of fl6| . These 
results extend the analytical understanding of the ana- 
log Hawking radiation to the sharp interface limit oppo- 
site to the hydrodynamic regime previously considered in 
the literature. Even if the singularity at the step makes 
the analogy with gravitational physics strictu sensu no 
longer valid, still the remarkable properties of the re- 
sulting emission support our choice of calling it analog 
Hawking radiation. 

The paper is organized as follows. In Sec. |TT]we present 
the physical system under consideration and we review 
the Bogoliubov description of quantum fluctuations in 
the presence of a sonic horizon. In Sec. IHII we derive ana- 
lytical formulas for the main observable quantities such as 
the emission spectrum and the density correlation func- 
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tion. These formulas are shown to successfully compare 
with the numerical results of [lq |. The effect of a finite 
initial temperature is discussed. Conclusions are drawn 
in Sec. [TV] 



II. THE PHYSICAL SYSTEM AND THE 
BOGOLIUBOV DESCRIPTION 

The system we consider is sketched in FigfT] A Bose- 
Einstein condensate is flowing along a one-dimensional 
atomic waveguide. For simplicity, the transverse trap- 
ping is assumed to be tight enough for a one-dimensional 
description to be accurate. The flow velocity is directed 
along the positive- a; direction, v(x) > 0. The density pro- 
file n(x) and/or velocity field v(x) of the condensate show 
a significant spatial modulation in the region around 
x = 0. Far from this interface region, both the density 
and the flow velocity tend to their up-stream (x — > — oo) 
and down-stream (x — > +oo) asymptotic values n u .d and 
v Ut d- The flow is kept stationary in time by means of a 
suitable external potential V ex t(x) and/or spatial mod- 
ulation of the atom-atom interaction constant g(x), so 
that the condensate wavefunction ipo(x) is a solution of 
the stationary Gross-Pitaevskii equation 
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V ex t(x) +g(x) | Vol 



(1) 

at a chemical potential fi. A dumb-hole configuration is 
realized when the up-stream region is subsonic v u < c u 
and the down-stream region is instead super-sonic Vd > 
Cd- Here, c u ,d — y guMnu,d/rn is the speed of sound and 
.9m, d = g{x — > — oo, +oo) are the asymptotic values of 
the interaction constant in the u, d-regions, respectively. 
For later purposes, we also introduce the corresponding 
upstream (downstream) healing length ^ n c ; = fr/(mc Ui d), 
chemical potential [i u ^d = mc^ d , and condensate wave- 
vector k Utd = mv u ,d/ 'fr- 
hi this dumb-hole configuration a sound wave emitted 
in the down-stream region is dragged away by the flow 
without being able to reach the up-stream region. The 
down-stream region is thus the sonic analog of the interior 
of a gravitational black hole, and the transition point 
where the velocity of the flow is exactly equal to the 
speed of sound is the analog of the event horizon. We 
conventionally locate this point at x = 0. 

Provided the condensate is everywhere dilute 
n(x)^(x) 3> 1, small fluctuations on top of the conden- 
sate can be described by the Bogoliubov theory of dilute 
condensates [l9j]. In particular, the elementary excita- 
tions correspond to the eigenvectors of the Bogoliubov 
operator: 
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A. Bogoliubov modes in a homogeneous system 

Within each homogeneous u, d region far from the in- 
terface, the eigenvectors of operator @ at an eigenenergy 
frui are plane waves of the form 



u(x) 
w(x) 



u(x) e lku ' dX 



(3) 

The wavevector k and frequency u satisfy the Bogoliubov 
dispersion relation u = Wg(fc) with 



\{k) = v u ,dk±c Uid k\ll + -(k£ Utd ) 



(4) 



As usual in Bogoliubov theory [19| , the upper (lower) sign 
in Eq. (j4|) refers to the positive (negative) norm branch, 
for which \U k \ 2 - \W k \ 2 = +1 (-1). 

Let us restrict our attention to the positive frequency 
modes (lu > 0) . The negative frequency modes can in fact 
be obtained from the positive frequency ones by simply 
reverting the sign of k and exchanging the values of the 
Bogoliubov coefficients U and W; corresponding modes 
then have opposite norms. 

In the up-stream sub-sonic region, for any ui > there 
are two positive norm modes satisfying the Bogoliubov 
equation lu = Wg(fc) with real wavevectors fcjf, A:° ut re- 
spectively. These modes correspond to propagating plane 
waves. In the Bogoliubov dispersion shown in Fig[2] 
(left panel), they correspond to the branches indicated 
as u-in and u-out. Throughout the whole paper, modes 
will be labelled as in-going ("in") if their group veloc- 
ity v g — diu/dk points toward the horizon and out-going 
("out") in the opposite case. In addition to the positive 
norm modes, a second pair of negative norm of modes ex- 
ist with complex wavevectors satisfying the Bogoliubov 
equation lu = w B (/;): these correspond to exponentially 
growing or decreasing evanescent waves [2(| ■ 

In the down-stream super-sonic region, two positive- 
norm real solutions of the Bogoliubov dispersion at k d \ 
and fc^i* exist for any ui > 0: In the Bogoliubov disper- 
sion of Fig[2] (right panel), these wavevectors correspond 
to the dl-'m and (il-out branches. On the other hand, 
two cases have to be distinguished in the negative norm 
sector, depending on whether lu is larger or smaller than 
the critical value f2 = w B (K) with K defined by 




(5) 



A pair of real solutions k l d n 2 and k°^ 1 to the negative-norm 
dispersion relation lu = lu^ (k) exist as long as lu < CI 
and correspond to the c?2-in and d2-out branches of Fig[2] 
(right panel). The critical frequency Q is the maximum 
of the d2-in and <i2-out branches. For lu > fl, these 
real solutions are again replaced by a pair of complex 
solutions corresponding to evanescent waves. Of course, 
the d2 branches do not exist if the condensate flow is 
everywhere sub-sonic. 
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FIG. 1: (color online) Sketch of the one dimensional dumb hole configuration considered in this paper. A BEC is flowing in 
the positive x direction. The horizon is located at x = 0. The flow is uniform in the asymptotic regions far from the horizon, 
with a velocity which is sub-sonic in the upstream region and supersonic in the downstream one. The wiggly arrows illustrate 
one among the different scattering processes described by Eq. (|10|l : an incident u-in wavepacket is partially reflected onto the 
u-out branch and partially transmitted on the dl-out and d2-out branches. The chosen labelling of the modes is explained in 
the text and illustrated in Fig0. 




FIG. 2: (color online) The Bogoliubov mode dispersion relations Q in the asymptotic upstream (left panel) and in the 
down-stream region (right panel). Solid (dashed) lines correspond to modes with a positive (negative) Bogoliubov norm. 



B. The scattering solution 

Far from the interface, the generic eigenfunction of the 
Bogoliubov operator {2]| at a frequency lu is built by su- 
perimposing within each u, d region all available propa- 
gating (that is non-evanescent) plane waves at the given 
frequency u>: 




(6) 



Here, plane waves have been separated into in-going 
(in) and out-going (out) ones according to the sign of 
their group velocity. Most remarkable among these so- 
lutions are the so-called scattering solutions, that de- 
scribe a plane-wave excitation originating from infinity 
(either upstream or downstream) on a well defined in- 
going mode, impinging on the horizon, and then leaving 
again towards infinity as a superposition of the different 
out-going branches: only a single in-going (3\ n amplitude 



is then non-vanishing, while two (if lu > 51) or three (if 
uj < fi) out-going /3° ut components have a finite value, 
describing reflected and transmitted waves. 

As a specific example, in the case where the in-going 
i = dl,d2 wave originates from the down-stream region, 
the scattering solution has the form 



u(x) 
w{x) 




(7) 



in the down-stream region and the form 



u(x) 



u kT , 



47T|t£S| 



(8) 



in the up-stream region. The two out-going components 
in Eq. (|7J) correspond to the reflected part in respectively 
the dl and d2 channels; of course, a d2 component can 
only be present if u> < tt. On the other hand, only a 
single outgoing transmitted component is present in the 
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upstream region. Similar expressions can be straightfor- 
wardly written in the case where the in-going u-m wave 
originates from the up-stream region. 

The expressions (0 and (JSj) based on the plane-wave 
expansion are only valid in the asymptotic regions far 
from the interface. For any given realization of the dumb- 
hole configuration the complete scattering solution for all 
x can be obtained by solving the full Bogoliubov problem 
@. A numerical solution for the case of a smooth inter- 
face was presented in [3. An analytical solution for a 
simplest step-like case will be presented in Sec lII CI 

In general, the linear relation between the amplitudes 
/3" 1 and /3° ut of in- and out-going waves can be written 
in a compact matricial form [7[. In the case of a generic 
flow that remains sub-sonic in both asymptotic regions, 
as well as for u> > fi in a dumb-hole configuration, this 
relation has the two-by-two form: 



= S(w) 



(9) 



which only involves positive-norm modes. On the other 
hand, the negative-norm d2 mode appears as soon as fre- 
quencies oj < Q are considered. In this case, the matricial 
relation has the three-by-three form: 



03? M 



= S(w) ( 0&(w) 

$&(<») 



(10) 



With the chosen normalization of the solution (jfjj), the 
S matrix connecting the (3 coefficients in the asymptotic 
regions is unitary in the Bogoliubov metric rj inherited 
by the norm of the corresponding plane-wave modes 



S(w) t 7/S(cj) = 



n 



(ii) 



More specifically, 77 = diag(l, 1,-1) in the case ui < ft, 
whereas 77 is just the standard 2x2 identity matrix r] = 
diag(l, 1) for u > O. 

Physically, the square moduli |Sj_,-| 2 of the S-matrix 
elements give the trasmission/reflection coefficients for a 
j-ingoing mode which scatters into an z-outgoing mode. 
The property [fTTjl ensures total energy conservation. It 
is interesting to note that a closely related approach was 
used in the context of quantum evaporation from super- 
fluid 4 He Hl|. 

As an illustrative example, for each quantum that in- 
cides on the system from the dl in-going mode with 
an energy htv, one has |S u( ii| 2 transmitted quanta in 
the u-out mode, | S^idi | 2 reflected quanta in the dl- 
out mode, and |S^2di| 2 quanta reflected in the d2-out 
mode. While the energy of the transmitted it-out and re- 
flected dl-out quanta is positive and equal to Huj, the 
energy of a quantum of the d2-out Bogoliubov mode 
is negative and equal to — fku 19]. As a result, en- 
ergy conservation then recovers the unitarity condition 

1 = S„dl + Sdidi — Srf2dl • 



C. The step like geometry 

All the theory reviewed so far is not limited to a par- 
ticular realization of dumb-hole and can be applied to 
generic configurations as long as the density and flow ve- 
locity tend to constant values in the asymptotic regions 
far from the interface. A numeric solution for a smooth 
interface was indeed reported in [l8[ . 

In the present paper, we restrict our attention to a 
model step-like structure with V ex t(x) = V^ t Q(—x) + 
V^ [t Q(x) and g(x) = g u Q(—x) + gd®(x) for which ana- 
lytical formulas for the components of the S-matrix can 
be found in closed form. As usual, 6 is here the Heavy- 
side step-function. Provided one imposes fl6| 



(12) 



the time-independent Gross-Pitaevskii equation (TT|) has 
a solution in the plane-wave form 



ipo(x) = y/nQ exp (ik x) 



(13) 



that describes a flow with uniform density n Ut d = no 
and velocity v Ui d = vq = hko/m. If one chooses Cd < 
vq < Cu, a dumb- hole configuration is obtained. The 
value of the potential jump ]/" t — V^ xt is then fixed by 
the condition (fT2"]) . This model configuration has the 
remarkable advantage of allowing for analytical insight 
without having to solve any differential equation. In a 
future publication we will present the analysis of other 
realistic dumb-hole configurations [22j]. 

In the step like geometry, the expressions ([7]) and ([8]) 
for the scattering solution can be extended to the whole 
x > and x < regions simply by including evanescent 
modes as well [2(| ■ A complete scattering solution is then 
formed by the superposition of five two dimensional col- 
umn vectors of the form ([3]): one of these corresponds 
to the incoming (i = u,d\, d2) mode and the four others 
are outgoing waves (two reflected and two transmitted). 
Depending on the value of w, some of the outgoing waves 
may be evanescent: in this case, between the two con- 
jugate complex wavevectors that solve the Bogoliubov 
dispersion ^ , only the wavevector value giving an expo- 
nentially decaying wave at large distance from the hori- 
zon has to be considered. 

Across the step-structure at x = the column vec- 
tor eigenfunction of C and its first derivative have to be 
continous. This provides the four matching conditions 
that are necessary to determine the amplitudes of the 
four outgoing waves as a function of the amplitude /3™ of 
the in-going i one. The coefficients of the non evanescent 
outgoing waves correspond to the matrix elements Sji 
with j — u, dl, d2. As usual, the d2 element only exists 
for us > il. Note that correct inclusion of the evanescent 
waves is here crucial to be able to fulfill the four matching 
conditions. 
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FIG. 3: Transmission and reflection coefficients for a dl in- 
going mode on a step-like dumb hole configuration. |S. u[ ji| 2 
(dotted curve) corresponds to the transmission on the it-out 
mode; | Sdidi | 2 (solid curve) and |S<j2di| 2 (dashed curve) cor- 
respond to the reflection on the dl-out and d2-out modes 
respectively. The parameters vp/cu = 0.75, vo/cd = 1.5 used 
here are the same as in Ref. [lfj]- Above f2 (Q/n u = 0.316 
for the chosen set of parameters) one recovers (insets a.l) and 
a.2) the usual case of a single transmitted (u-out) and a sin- 
gle reflected (dl-out) waves, respectevely. Inset b) displays a 
comparison of the exact |Sdidi| 2 coefficient (solid line) with 
its low energy approximation (TH)) (dot-dashed line) . 



D. Transmission and reflection 

The typical behaviour of the scattering coefficients for 
a dl ingoing wave is plotted in Fig{3] as a function of the 
frequency lo. For ui > il, the transmitted u-out wave and 
the reflected (il-out wave are only present. As usually 
expected in wave mechanics, the transmission (reflection) 
coefficient increases (decreases) as a function of u>, see the 
insets (al-a2). Energy conservation imposes that the sum 
of these two coefficients is always equal to unity. 

For uj < fl also the (negative-norm) reflected mode 
d2-out is involved in the dynamics. Remarkably, all 
three scattering coefficients diverge as l/u> in the low 
lu limit. Nonetheless, energy conservation is ensured 
by the 77-unitarity of the S-matrix, which now imposes 

I'S'udll 2 + ISdldll 2 — |<Sd2(il| 2 = 1- 

Straightforward algebraic manipulations lead to a sim- 
ple analytical formula for the reflection coefficient in the 
low-w limit. 



>dldl 



c« (Cu - vq)(c u - c d ) 
Cd {c u + v )(c u + c d ) 



3/2 



(14) 

the excellent agreement between this low energy approx- 
imation and the exact Bogoliubov result is apparent in 

Fig. EP>). 

The wave mechanics encoded in the S-matrix is physi- 
cally illustrated in FigH]where we show the result of a nu- 
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FIG. 4: A dl ingoing wave packet is incident on a step-like 
sonic horizon located at x — (upper panel). When the 
central frequency omega of the wave packet is smaller that 
Omega, a single transmitted u wavepacket and two dl, d2 re- 
flected wave packets are visible (middle panel). A single u 
trasmitted and a single dl reflected wave packet are instead 
visible for u > Q. (lower panel). Same system parameters as 
in FigO Wave-packet carrier wave-vector q^ u — —1.2 (upper 
and middle panels), q£ u = —1.35 (lower panel). 



merical solution of the time-dependent Gross-Pitaevskii 
equation that describes the scattering of a incident dl 
in-going wave packet (upper panel) onto the black hole 
horizon. For a low initial energy, the in-going wave packet 
splits into a u-out transmitted one (central panel), and 
a pair of dl, d2 reflected ones: Among these, the bigger 
and slower wavepacket corresponds to the negative-norm 
d2-out mode, while the other one corresponds to the dl- 
out mode. As expected, for an in-going wavepacket en- 
ergy above £1 (lower (c) panel), only the u transmitted 
and the dl reflected wavepackets are visible. A related 
phenomenology was experimentally observed for surface 
waves in a water tank in [lOj and theoretically discussed 
in [l3| for a chain of trapped ions. In j24[, this nega- 
tive norm reflected mode was interpreted in terms of a 
bosonic analog of Andreev reflection. 



E. Quantization of the modes 

In the presence of a macroscopically occupied conden- 
sate, the full Bose field operator can be expanded as 



%j}{x) = ipo(x) + S'tp(x) . 



(15) 



The coherent condensate is described by a classical field 
ipo that evolves according to the Gross-Pitaevskii equa- 
tion. Quantum fluctuations are described by the opera- 
tor term Stp whose dynamics is ruled by the Bogoliubov 
operator 

A most favourable basis to our purposes consists of the 
scattering solutions at a given frequency u> that we have 
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discussed in Sec lIIBI Within each uj sub-space, the three 
scattering modes corresponding to the / = u 1 dl 1 d2 in- 
going modes form an ortho-normal basis (of course, one 
has to restrict to I = u, dl for uj > Q) in the Bogoliubov 
fy-metric, 



dx (u Iul (x)*u Ju! '(x) -w} u (x)wj u >(x)) 

= siSjjS(uj — uj') (16) 
dx {u Iul (x)v Ju; >(x) - v Iu] (x)u JuJ >{x)) = . (17) 



The si coefhcient gives the sign of the Bogoliubov norm 
of the mode: in our case, one has sj = +1 for I = u,dl 
and sj = — 1 for / = dl. Spatial integration is over the 
whole space. 

The final form of the quantum field is obtaine d by 
the standard Bogoliubov quantization prescription [l9j: 
the amplitudes corresponding to positive-norm modes be- 
come destruction operators while amplitudes of negative- 
norm modes become creation operator, 

/>oo 

S'i/j(x) = I du) 1 ^ [uiuj(x) aj(u) + w* Iu (x) a/(w) t ] 
Jo T _„ 



I=u,dl 



\u d2 ^{x) d\ 2 {uj) + w* d2 u) {x)d d2 {ui) . (18) 



Expectation values of physical observables (e.g. the den- 
sity correlation function) are then straightforwardly eval- 
uated by imposing the suitable boundary conditions on 
the expectation values of products of in-going operators, 
di and a j. 

If one is interested in correlation functions involving 
the out-going modes only, a reformulation of (fTS)) in 
the input-output language [23[ can be used, as proposed 
in [3]. In our case, the input-output relations consist of a 
linear relation connecting the operators of the out-going 
modes b Ut di,d2 to the in-going d u ^di,d2 ones via the S ma- 
trix (|10[) introduced in the previous subsection: 



b u (ui) 
bdi(w) 



SM 



a u (uj) 



(19) 



As a consequence of their negative Bogoliubov norm, the 
d2 modes appear in both the left- and the right-hand side 
of (|19p as creation operators rather than destruction ones. 
As we shall see in Sec lIII Al this simple fact is the mathe- 
matical origin of the Hawking emission in our formalism. 
From a quantum optical perspective, Hawking emission 
can then be interpreted as parametric down-conversion 
of Bogoliubov sound waves by the horizon. 



III. OBSERVABLES 

The expression (IT51) for the quantum field in the Bo- 
goliubov approximation and the input-output relation 



([19]) are the starting point for the calculation of physical 
observables of the system: in particular, our attention 
will be focussed on two most remarkable ones, namely 
the spectral distribution of the Hawking emission and 
the long-distance behaviour of the correlation function 
of density fluctuations. 



A. Emission spectrum 

In the gravitational context, the only observable quan- 
tity is the Hawking radiation outside the black hole. One 
of the most remarkable feature is that this radiation is 
thermal at a temperature univocally determined by the 
surface gravity of the black hole. 

In our condensed matter context, the Hawking emis- 
sion from the dumb hole corresponds to the phonons that 
are emitted into the upstream region on the it-out branch. 
Assuming that no correlation between the different in- 
going modes exist, the emission spectrum, that is the 
number of phonons emitted per unit time and per unit 
bandwidth into the u-out branch is straightforwardly cal- 
culated from the input-output relation (|19p : 



dl?. 



dt duj 



= (bi(u)b u (oj)) = 



|Swi| 2 Idi 



>ud2\ 



(I 



(12 



1) 



(20) 



For a system initially at zero-temperature, all the 
(to) a Ut di,d 2 {^>)) = 0, yet a finite Hawk- 



,,dl,d2 



l u.dl,d2 



ing emission exists as a consequence of the +1 term aris- 
ing from the (d d 2(uj) a d2 (uj)) expectation value that en- 
codes quantum fluctuations. 

Restricting our attention to the low-frequency part of 
the spectrum, an analytical form can be found for the 
emission spectrum: 



dl?, 



dt duj 



>t(l2 



HI 5 



T=0 

(c« 



v ) 



■ V ) 



3/2 2u 



It is remarkable to note that this formula still gives the 
typical the 1/lu thermal behavior of Hawking radiation 
even though one is not allowed to use the gravitational 
analogy: in the present step-like case, the surface gravity 
is in fact formally infinite. The effective temperature of 
the emission, i.e. the coefficient of the l/uj term, is here 
determined by the microscopic physics of the condensate 
that fixes the cut-off frequency Q and the chemical po- 
tential [l u ,d- 



B. Density correlations 

As first remarked in Ref.[l5|], the density correlation 
function appears to be the most promising tool for iden- 
tifying Hawking radiation: in particular, this quantity 
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was at the heart of the numerical observation of [l6| . An 
example of such a calculation is reproduced in Figj5] The 
dark regions correspond to antibunching, and the bright 
ones to bunching. The strongest feature in this figure is 
the dark stripe along the x — x' line: as its origin is well 
known and corresponds to the antibunching due to the 
repulsive inter-atomic interaction [25j , we will not discus 
it further. In the present paper, we will rather focus our 
attention onto the other features of the correlation plot, 
labeled as u — d2, u — dl and dl — d2, that encode the 
information on the Hawking emission. 

Within Bogoliubov theory, the density correlation 
function in the stationary state 



G^(x,x') = 



(ft(x)ft(x')ip(x')4>(x)) 



{^{x)4i{x))(^{x')^{x')) 
can be expanded in its w components as: 



- 1 



G^{x,x') 



G i2) (u,x,x')duj 



(22) 



(23) 



Each u> component has to be evaluated imposing suit- 
able boundary conditions on the expectation values of 
in-going operator averages. In the following of the discus- 
sion, it will reveal useful to separate the zero temperature 



(2) 

contribution Gq from the thermal one G th 

G (2) (lu, x, x 1 ) = (w, x, x') + G$ (lu, x, x') . (24) 

The former term only involves the zero-point fluctuations 
of the in-going modes 

n G^\uj,x,x') = — [ w* Iu ,{x)r IuJ {x) 

I=u,dl 

+ U*d2u>( x ) r d2u J (x')+C.C. , (25) 

while the latter includes the initial populations 
(a\(uj) ai(uj)): 

noG < $[u,x,x') = — Y [r* Iu (x)r Iu (x') + c.c] x 



(2) 



.dl,d2 



X (4(cj) ai(u)) . (26) 
In both these formulas, we have set rj = uj + wj. 

C. Zero Temperature density correlations 

1. In-Out correlation 

Let us first consider the in-out zero-temperature cor- 
relations between the inside (x' > 0) and the outside 
(x < 0) regions. For x and x 1 far from the interface, 
it is legitimate to only keep the terms with a stationary 
phase and neglect the fast oscillating ones. Making use 




FIG. 5: (color online) Plot of the rescaled density correlation 
no £ u [Gq 2 ' (x, x') — 1] for an initial zero t emp erature. The 
data are the same as the one presented in [lq |. but the color 
scale has been modified in order to make the positive u-dl 
correlation clearly visible. This signal was overlooked in the 
original paper given its relative weakness. The dashed lines 
identify the different feature of the correlation function as 
discussed in the text. The white dotted line corresponds to 
the cut x' — x — 39.5 £ u used in Figs. [7] and [6] 



of the r^-unitarity of the scattering matrix and setting 
Rk = Uk + Wk , the simple expression 



n fl G{j 2) (w,i,3;') = ^ [ 



iJ/cout (w)i? fe out (u>) 



l=d\.d2 



,,out ,,out I 
u g,u VgA I 



e *K u V)*-fcrM*'] S : d2 M S ld2 (cu) + c.c.}, (27) 

is obtained for the long-distance zero-temperature corre- 
lation function. As all the S-matrix elements appearing 
in Eq. (|27p involve the d2 mode, it is apparent that a 
non-zero long-distance correlation is only possible in a 
dumb-hole configuration: in terms of the input-output 
picture of (|19[) . the two terms in ([27)) originate from the 
quantum fluctuations of the ingoing d2 mode that are 
scattered in either the u or the I = dl, d2 modes. 

As the integrals in (|23p are dominated by the small oj 
part of the integral, simple analytical expressions can be 
worked out by extrapolating the low-w asymptotic form 
of the S-matrix elements 



K./2 



M Si d2 (u)) 



(28) 



to the whole spectrum and then cutting the integral at CI. 
In this hydrodynamic approximation, fc° ut ~ ^/Vj (j = 
u, dl, d2) and the velocities Vj of the outgoing modes are 
V u = vq - c u < 0, V d i =v + c d > 0, V d 2 = v - c d > 0. 
As typical of Bogoliubov theory [19| , the Rk coefficients 
giving the density weight of the Bogoliubov modes scale 
as \fuj. The in-out correlation function then takes the 



simple form 



n G^\x,x') 



l=dl.d2 



A l cosine ~ j)) 



(29) 



where sine (2;) = sinx/x and the Ai coefficient has the 
following explicit form 



i-dl (d2) 




3/2 



Cd + (-)c ll 



in terms of the microscopic parameters of the system. 

From the explicit expression (|29|) . it is immediate to 
see that the density correlation is maximum along the 
pair (I = dl, d2) of straight lines 



x x 

Vu~ = vi 



(31) 



that exit in opposite directions x < 0, x' > from the 
horizon position x = x 1 = 0. The absence of lateral shift 
from the horizon position x = x' = is a consequence of 
the reality of the Ai coefficients. As a consequence of the 
inequality c u > vq > Cd, the sign of the d2 contribution is 
positive, while the one of the dl contribution is negative. 
The d2 contribution is always stronger than the dl one 
by a significant factor 



\A d 2/Vd2\ _ (C u - C d ){vo ~ C d ) 

\A d i/Vdi\ (cd + Cu)(vo + Cd) 



(32) 



The d2 contribution corresponds to the Hawking signa- 
ture anticipated in [l5[ and numerically observed in . 
Once the horizon is formed, the quantum fluctuation of 
the incoming negative-norm d2 mode start being con- 
verted into correlated pairs that emerge from the horizon 
in the u and d2 modes. The two phonons are simulta- 
neously created at the horizon and then propagate away 
at speeds V u — vq — c u < and Vd2 = vo — c d > 0. At 
time t after their emission, they are therefore located at 
% = V u t < and x' = Vd2t > 0, which explains the corre- 
lation between the density fluctuations at points verifying 
x/V u =x'/V d2 . 

This simple argument can be straightforwardly ex- 
tended to explain the geometry of the other dl feature. 
This feature went overlooked in |l3 because of its weak 
intensity, but was mentioned in [181 ]. When the numeri- 
cal data of [lj| are plotted with a suitable color scale as 
done in FigO the u-dl correlation is perfectly visible as 
a positive correlation tongue. 

A quantitative comparison of the prediction (|27p with 
the numerical results is illustrated in Fig|6] where cuts 
of the density correlation function along the line x' — 
x = 39.5 £ u (marked by a white dashed line in Fig[S]) are 
shown using the same set of parameters as in Ref.[lj|. 
In spite of the approximations made, the agreement is 
remarkable. 

The accurateness of hydrodynamical approximation 
(|2"9")1 is finally validated in Fig[7] As shown in the main 




20 40 

(x+x')/2^ 



FIG. 6: Cut of the zero-temperature density correlation func- 
tion Gq alone the x' — x — 39.5 ^ u line. Solid line: numerical 
result from [lq . Dashed lines: prediction of Eqs. (|27p and (|33p 
for respectively the u — d2 and the dl — d2 contributions. On 
this scale the u — dl contribution is hardly visible. 



a 



0.004 




1 1 


1 1 1 


1 1 


0.002 





















-0.002 










-0.004 












"0.0006 
" 0.0004 


1 1 1 1 1 1 1 

- a) A- 


\ I / / 


1 1 1 /Is. 1 1 1 




-0.006 


"0.0002 





\ \ 1 -0.0005 
\ A / -°' 001 






-0.008 


"0.0002 







\ \ l I 1 1 1 ~ 






-10 ! 10 2 


20 3p 40 50 60 







FIG 

f(2) 



Cq' along the x' 



(Color online) Cut of the density correlation function 
x = 39.5 £ u line. Black solid lines are 
the prediction of (|27[) and (|33[) . Red dashed lines are the hy- 
drodynamical approximations (|29[) and (|36[l . The main panel 
shows the u — d2 feature in the in-out region. Inset (a) shows 
the u — dl feature in the in-out region. Inset (b) shows the 
dl — d2 feature in the in-in region. 



panel, the agreement for the I = d2 feature is quite 
good even though the hydrodynamic approximation is 
not fully able to reproduce the asymmetry of the peak 
nor its slightly back-shifted position with respect to the 
straight line (|3ip . Both these facts can be traced back to 
the quick deviation from a linear dispersion that is visi- 
ble in Figl^l (right panel) for the d2 Bogoliubov branch. 
This interpretation is confirmed by the excellent agree- 
ment which is instead found for the I = dl contribution 
(inset (a) of FigJTJ) : the discrepancy of the dispersions of 
both the u and the dl modes from the hydrodynamic ap- 
proximation is in fact very small in the whole frequency 
range up to w = fi. 
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2. In-In correlation 

The in- in correlation function for a pair of points x, x' 
both located inside x, x' > the dumb hole can be calcu- 
lated along the same lines. In this case, a single contribu- 
tion appears that originates from the correlation between 
the dl and the d2 modes, 



n G 2) (w, x,x') 



1 ( RwMRk&iu) 



4tt 



I,, out „,out I 
\ U n,dl U a,d2\ 



J g,dl u g,d2\ 

x e <(*3S t M«-fc3S t (<-K) S* dld2 (cj) S d2d2 (co) + c.c. 

+ (x <-> a/). (33) 

Even thou gh it did not appear in the analytical calcu- 
lations of 15( based on the gravitational analogy, this 
contribution was observed and discussed in the numeri- 
cal paper 

The same hydrodynamical approximation that led to 
(|29|l can be used to obtain an approximate expression to 
this dl — d2 feature as well. For small w, the S-matrix 
product can be expanded as 



BjJ, u 
fko 



with 



B 



Cu (C u ~ Vq) C d \ 



2c d (c„ + v ) \cl 6, 



(34) 



(35) 



Extrapolating this expression to the whole u> spectrum 
and cutting the integral at Q, one immediately gets to 
the expression 



n Q G^\x,x') 



Bncl sinc(rj(# 



+ {x^ x') . (36) 



D. Finite Temperature density correlations 



vq and sound velocity c u and is at thermal equilibrium in 
the moving frame at a temperature T [261 ]. The horizon 
is then (adiabatically) switched on by ramping down the 
scattering length in the downstream region. 

Within each semi-infinite uniform region, the popu- 
lation of each k mode is adiabatically preserved dur- 
ing the formation of the dumb hole. The initial popu- 
lation (ajaj) of each / = u,dl,d2 mode at wave vec- 
tor ki is given by the thermal population u(Oj) = 
[exp(Mli /HbT) — corresponding to its frequency 



h 2 kj (h 2 k 2 T n „ 



(37) 



in the comoving frame before the creation of the horizon. 
For each ui, the wavevector k d of the corresponding in- 
going mode has to be determined by inverting the relation 

The density correlation function at a finite tempera- 
ture can then be written as a sum of different terms, 

G (2) (w, x, x 1 ) = G 2) (co, x, x') [1 + n(n dl ) + n(fl d2 )] 

+ F sh (oj,x,x') + F nsh (uj,x,x'). (38) 

Here, n(u>) = [exp(hui/kBT) — 1] _1 is the thermal law at 
a temperature T . 

The first term in the r.h.s. of (|38|) is straightforwardly 
interpreted as a thermal enhancement of the zero-point 
signal described in the previous subsection, i.e. a stim- 
ulated Hawking emission. Other scattering processes pe- 
culiar to the T ^ case are responsible for the other 
contributions, which can be separated in a u — d2 term 
due to the presence of a sonic hole, 



n F S H(uj,x,x') 



„OUt,,.OUt I 

g,u u g,d2\ 



S* uu Sd2u(n(n u ) - n(n dl )) + c.c. , (39) 

and a u — dl term that instead persists even in absence 
of the sonic hole: 



One of the major issues in view of an experimental 
observation of the analog Hawking radiation is the role 
of temperature. To this purpose, a naive, but extremely 
constraining requirement is often imposed that the tem- 
perature of the system should be (much) smaller than 
the Hawking temperature of the emission. The numerical 
simulations of 16] have instead shown that correlations 
are robust with respect to temperature, which supports 
once more their promise in view of an actual experiment. 
In the present section we extend the theory of correla- 
tions developed in the previous section to the case of a 
condensate at a finite initial temperature. 

In particular, we focus our attention on the same dy- 
namical situation already considered in 16]: The con- 
densate initially has a uniform density no, flow velocity 



n F NS H(u,x,x') 



i? fc out R k a 



i(fe° 



471-4/ k 



J g,u g,dl I 

S;„ S dlu (n(Q u ) ~ n(O d i)) + c.c. . (40) 

In the low-w limit the S-matrix element products in- 
volved in the new terms (l39l) and (l40l tend to finite values 




c d + (-)t'o 



Cd c u + v c u + v 
while the thermal population is proportional to 1 /lo 

Vq + Cu k B T 



(41) 



n(f2„) - n{Q dl ) 



(42) 
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Combining these two facts, it is immediate to see that 
in the hydrodynamic approximation the contributions of 
(|39|) and ([40]) to respectively the u—dl and u—d2 features 
then the same form as the zero-temperature one (I29[) . 

Summing up all the terms, the only effect of a finite 
initial temperature on the dl — d2 feature is a bosonic 
stimulation factor 1 + n(Q<nJ + n{p.d2)- In particular, 
this feature is not affected by the Fsh and Fnsh terms. 

The u — d2 feature is affected by this multiplicative 
factor as well as by the additional term Fsh'- remark- 
ably, the zero-point (f25)) and thermal ([55)) contributions 
have opposite signs. Within the hydrodynamic approxi- 
mation, the ratio of their contributions to the peak cor- 
relation signal has the simple expression 

_ c u (v - Cd){c u - c d ) k B T 

Oo - c d) 3/2 Mu ■ [6) 

As usual, the zero-point contribution dominates at low 
temperature while it is overwhelmed by the thermal one 
at high temperature. For the parameters used for the 
finite temperature numerical simulations in Ref . [l6| , the 
T = was still the most relevant one. 

The effect of the thermal occupation on the u — dl fea- 
ture is even more dramatic. Also in this case, the (am- 
plified) zero-point positive feature due to the I = dl term 
in (|29|) and the thermal contribution (|40|) have opposite 
signs. As the involved scattering process is not limited to 
the uj < ft window, the term (|40p gets contributions from 
high frequencies and is for this reason generally larger in 
magnitude than (|39|) . For this reason, the zero-point con- 
tribution only dominates at very low temperatures and 
is quickly overwhelmed by the thermal one: this effect 
is clearly visible when comparing the T = numerical 
simulations shown in FigJH] with the finite temperature 
ones shown in Fig. 6(a) of Ref. [To] ]. 

When the flow remains everywhere subsonic, all long- 
distance features of the density correlation function dis- 
appear but for the positive u — dl tongue at finite tem- 
perature. This fact was apparent in Figs. 4(a) and 6(b) of 
Ref. [l6| and corresponds within the present formalism to 
the Fnsh term. Its physical origin can be traced back to 
the finite reflectivity of the interface region even in the 
absence of an horizon and the different initial occupation 



of the u and dl incident modes. 

Before concluding, it is important to stress that a dif- 
ferent configuration was considered in [Lsl ]. where the 
dumb hole was assumed to be in thermal equilibrium 
in the comoving frame after completion of the horizon 
formation process. Of course, this configuration can be 
still described by Eq. ((38|) . but the initial populations 
n{£ldi,d2) have to be computed using in (|37p the Cd in- 
stead of c„ for the downstream region. 



IV. CONCLUSION 

In this work we have made use of the Bogoliubov the- 
ory of dilute Bose-Einstein condensates to study the ana- 
log Hawking emission that is emitted by an acoustic black 
holes. Our framework does not rely on any gravitational 
analogy and is able to provide closed-form analytical for- 
mulas for the emission spectrum and the density corre- 
lation function of a simplest step-like configuration. Al- 
though the surface gravity is in this case formally infinite, 
the low-energy part of the emission spectrum still shows 
a thermal character, but the temperature is found to de- 
pend on the microscopic properties of the condensate. 
The extension of our theory to the case of a non- vanishing 
initial temperature is discussed. We have shown that this 
general framework is able to quantitatively reproduce the 
results of first principle numerical calculations [16j . and 
to provide a physical understanding of the various fea- 
tures that appear in the density-density correlation func- 
tion. Generalization of the present theory to more gen- 
eral dumb-hole configurations is in progress [22| . 
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